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We compute primordial light-element abundances for cases with fine structure constant a different 
from the present value, including many sources of a dependence neglected in previous calculations. 
Specifically, we consider contributions arising from Coulomb barrier penetration, photon coupling 
to nuclear currents, and the electromagnetic components of nuclear masses. We find the primordial 
abundances to depend more weakly on a than previously estimated, by up to a factor of 2 in the case 
of ^\A. We discuss the constraints on variations in oi from the individual abundance measurements 
and the uncertainties affecting these constraints. While the present best measurements of primor- 
dial D/H, ■'He/H, and '^Li/H may be reconciled pairwise by adjusting a and the universal baryon 
density, no value of a allows all three to be accommodated simultaneously without consideration of 
systematic error. The combination of measured abundances with observations of acoustic peaks in 
the cosmic microwave background favors no change in a within the uncertainties. 

PACS numbers: 26.35.-|-c,24.10.-i,12.20.-m,98.80Ft,98.80.Es 



I. INTRODUCTION 

The recent literature contains claims for time variation 
in the fine structure constant a at a level of one part in 
10^ since a redshift of 2, based on observations of metal 
ions backht by distant quasars ||, |, || . Regardless of 
theoretical prejuduce for or against these claims ^, it 
is of interest to examine other constraints on a variation. 
In particular, changing the fine structure constant would 
influence the abundances of light nuclides produced dur- 
ing big-bang nucleosynthesis (BBN). BBN took place at 
redshifts 10^ < z < 10^°, ending minutes after the big 
bang, and it is the earliest cosmological test for which 
the basic physical laws are well-understood. The physics 
is less certain in varying-a scenarios, but BBN remains 
the earliest reasonable probe of variations in a. It thus 
provides a large lever arm in the sense that if a varies 
with time, larger variations may plausibly (though not 
necessarily 0, |8j) be observed over longer times. 

Past BBN calculations for varying-a models focused 
on the abundance of ^He In some ways, this 

keeps the physics simple, because a should affect the '*He 
mass fraction Yp only through the neutron/proton ratio 
at freezeout of weak interactions before BBN - any effect 
on the rates of subsequent nuclear reactions is irrelevant 
because essentially all neutrons are gathered into ''He 
nuclei by the end of BBN. The likely size of the effect 
on Yp is also large compared with observational errors 
if \5\ > 0.05, where S = {a — ao)/ao and subscript zero 
denotes the present value of a quantity. 

More recently, Bergstrom, Iguri and Rubinstein 
(hereafter BIR) extended the varying-alpha calculation 
to include other light nuclides by estimating the depen- 
dence of the charged-particle reaction rates on Coulomb 
barrier penetration and hence on a. However, they chose 



to neglect several places where a enters the rates (from 
both barrier penetration and other considerations), un- 
der the assumption that the exponential a dependence 
in the Gamow factor will dominate all other a depen- 
dence. Our work corrects these omissions, most impor- 
tantly those arising from barrier penetration and from 
coupling of photons to nuclear currents. We also esti- 
mate the effects of final-state barrier penetration and the 
a dependence of the nuclear masses. We use more elabo- 
rate models for the a dependences of the mass-7 radiative 
captures, and we correct BIR's use of reaction rates that 
were already out of date when they did their work. 

Finally, we elaborate the BIR discussion of BBN con- 
straints on a, examining the region of the (1^5/1^,0:) 
plane allowed by each of the observed light-nuclide abun- 
dances, (ris is the universal baryon density in units of 
the critical density; h is the Hubble constant in units 
of 100 km s~^/Mpc.) We then discuss the prospects 
for producing a consistent solution, given systematic un- 
certainties in the observed abundances. We leave aside 
the question of how (or whether) other parameters of 
the standard model of particle physics vary as a varies 
§, |ll|, |l|, |l§, though this may also affect BBN. 



II. LIGHT NUCLEI AND THE FINE 
STRUCTURE CONSTANT 

The BBN light-element abundances, apart from ^He, 
depend on the fine structure constant through nuclear 
reaction rates. Calculating how these rates depend on a 
is the main task undertaken by BIR. 

Charged particles must tunnel through a Coulomb bar- 
rier to react at low energies, and this produces most of 
the energy dependence in the cross section. Changing a 



2 



raises or lowers the Coulomb barrier, altering the reac- 
tion rates. To describe laboratory results and compute 
thermally-averaged rates, it is useful to separate out the 
Coulomb part of the energy dependence by writing the 
low-energy cross section as 



(1) 



Here, rj — aZiZ2y^ fJ,/2E, fx is the reduced mass, and Zi, 
Zi are the charges of the incoming nuclei. The astro- 
physical 5'-factor SijS) is then usually a slowly-varying 
function of energy. 

BIR calculate the a dependences of the thermal re- 
action rates that result if the S'-factors are assumed to 
be independent of a; i.e., if all of the a dependence is 
expressed in Eqn. |I and the definition of r\. We have ver- 
ified that their resu ts are correct, given this assumption. 
However, Eqn. |l| does not explicitly contain all of the 
expected a dependence, even if the only effect of chang- 
ing alpha is to change the Coulomb barrier. Moreover, 
BBN takes place at relatively high energies (a few hun- 
dred keV) and Z\Z2 is only in the range 1-4 in BBN, so 
the exponential does not dominate sufhciently to make a 
dependence of the S'-factor completely negligible by com- 
parison. If BBN involved energies a factor of two lower 
or nuclei with only slightly larger charges, corrections to 
the BIR treatment would be negligible. 

The first source of a dependence in S(E) is Coulomb 
barrier penetration. Cross sections are proportional to 
the barrier penetration factor vi: 
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Here, Fi and Gi are the regular and irregular Coulomb 
wave functions, respectively I is orbital angular mo- 
mentum of the entrance channel, k is the wavenumber 
in the entrance channel, and R is the nuclear interaction 
radius. (A clear derivation of Eqn. |2| may be found in 
Chapter VIII of Blatt and Weisskopf Q). Assuming 
Z = (a reasonable assumption at very low energies), v 
is nonzero even when i? = 0. For small fc, it is reasonable 
to approximate the barrier-penetration part of the cross 
section as 
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We have thus "derived" the heuristic expression in Ec^n. 
|l|, and found that S{E) is proportional to a. If the en- 
trance channels of all reactions are reasonably approxi- 
mated as S waves penetrating to zero separation, it fol- 
lows that the charged-particle rates in BIR should be 
multiplied by an additional factor of (1-1-5) = aBBN/c^o 
(where q:bbn is the value of a at the time of BBN). When 



R > OT I > 0, Eqn. ^ almost always provides an over- 
estimate of the sensitivity of Eqn. ^ to a. (Exceptions 
with ^ 10% underestimates of dvi/da occur &t E > 100 
keV for reactions with Z1Z2 — 1.) 

The pedagogical literature usually motivates Eqn. |l| 
by a WKB approximation. Even where WKB pro- 
vides a quantitatively good approximation (not the usual 
case [|l^ ) , the factor of a on the front of Eqn. ^ enters the 
WKB barrier penetration factor through the normaliza- 
tion that matches solutions across the classical turning 
point. 

When a reaction produces two charged particles, they 
must escape from a Coulomb barrier, and this introduces 
more a dependence in the rates. The cross section is 
again suppressed by a factor of the form shown in Eqn. 
^, but with E ^ Q + E and appropriate charges. The 
effect is weak because Q (the energy released) is gen- 
erally much larger than the Coulomb barrier, but there 
are two exceptions, '^He(n,p)'^H and ^Be(n,p)^Li. These 
two cases are very interesting, because their Q values de- 
pend strongly on a (see below), and this introduces a 
dependence in the final-state barrier penetration that al- 
most cancels the dependence due to the changing barrier 
height. 

The generally high Q-values in BBN invalidate the 
approximation in Eqn. ^ for final states. In fact, the 
final-state barrier penetration factors for the BBN reac- 
tions are almost independent of energy over the BBN 
energy range. Even Eqn. |^ badly overestimates the ef- 
fect because k is larger for the final states. We esti- 
mate the final-state suppression by assuming i? = 2 fm 
and / = 0,1,2, computing logarithmic derivatives of vi 
with respect to alpha, and looking for the largest effect. 
The results should be rough but reasonable upper limits 
to the size of the effect. To keep thermal rates simple, 
we exploit the weak energy dependence of these deriva- 
tives and multiply the rates by energy-independent fac- 
tors (I + aS) which model their a dependence. In these 
factors, a is estimated to be slightly larger in magnitude 
than avf^dvi/da for all I considered, and the adopted 
values are shown in Table |[ In the cases of '^He(n,p)^H 
and ''Be(n,p)''Li, the low Q-value leads us to include in 
a the complete dependence of final-state channel widths 
on Q. (Channel widths are proportional to the product 
of channel velocity and vi [|8) , both of which vary with 
Q.) This adjustment is based on our estimate of the elec- 
tromagnetic contribution to the masses (see below) and 
takes the corresponding numbers in Table | from slightly 
negative to slightly positive. 

The matrix elements for radiative captures - photon- 
emitting processes - are proportional to a because a is 
the strength of the coupling between the photon and 
nuclear currents. Therefore, all radiative capture rates 
should be multiplied by an additional factor of (1 -I- (5). 

Some of the radiative captures PH(q;, 7)^Li and 
■^He(a, 7)''Be] are also external reactions (e.g. |l^, pOf), 
so particles need not penetrate the Coulomb barrier all 
the way to the nuclear surface; the reaction matrix ele- 
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Reaction a Reaction ao ai as 04 05 

d{d,p) -0.16 ^H(a,7)'Li 1 1.3718 0.502 0.1829 0.2693 -0.2182 

3He(n,p) +0.30 ^He(a,7)^Be 1 2.148 0.6692 -5.566 -10.63 -5.730 

^He(d,p) -0.09 ^^^^^^^^^^^^=^^^^^^^^^^^^= 

^Li(p,Q) -0.18 

^Be(n p) +0 20 TABLE II: Estimated dependence of radiative capture rates 

^^^^^^^^^^^^^^^^^^=^^^^^^^^^^^^^^^^^^= on a, with /(q) — ai<5' as defined in Eqn. 0. 



TABLE I: Corrections for final-state interactions for processes 
involving two charged particles in the final state. To correct 
the rates, multiply by (1 + aS), where a is given above. 



merits contain large contributions from outside the nu- 
clear interaction radius. Further, Coulomb forces govern 
the forms of the final states (^Li and ^Be bound states) at 
these distances. Most importantly, the electrostatic en- 
ergies of the bound states are also large fractions of the 
reaction Q values, so that variations in a change the out- 
going photon energy significantly. This energy enters the 
cross sections a.s a E!^. We used a simple radiative cap- 
ture model to estimate the combined size of these three 
effects, though we found that the last one dominates. 

We assume the following: 1) These reactions can be 
treated as direct captures in which the bound states con- 
sist of smaller clusters: ''He and in the case of ""Li, or 
^He and '^He in the case of ^Be. 2) Wave functions can 
be generated from potential models; we use the potential 
of Buck et al. 3) The only effect of changing a is to 
change the charges in the Coulomb term of this potential. 
There is no question that the assumptions numbered 1 
and 2 are reasonable descriptions of laboratory data. It 
is plausible that assumption number 3 is no worse than 
other assumptions we have made implicitly in modelling 
BBN with varying a. We estimate the variation of the 
thermal reaction rates with a from this model. We find 
the effect on the thermal rates (beyond the effect in Eqn. 
|I|) to be roughly independent of temperature, and we 
characterize it by fits to the computed variation beyond 
the corrections already discussed. 
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(with (av) the thermally-averaged rate). These fits give 
estimates of the rates at a ^ ao that are accurate to 20% 
in the additional effect beyond the exponential factor and 
two powers of a already considered (at temperatures less 
than 10^" K), and they are shown in Table |l[ We assume 
that this a dependence holds relative to the reaction rates 
computed from laboratory measurements. (See Sec. Ill 
for the baseline a = rates.) 

Nuclear masses depend on a through electromagnetic 
interactions between nucleons. In turn, differences be- 
tween initial-state and final-state masses are used in nu- 
cleosynthesis calculations to determine the rates of re- 
verse reactions (by detailed balance). BIR report that 
reverse reactions have essentially no effect on the BBN 
yields, but we find that when we alter the masses, the 



TABLE III: Electromagnetic contributions AAIem to the re- 
action energy yields, from the calculations of Refs. |2^, 
Shown for comparison are the total energy yields Qo- 



Reaction AMem 


(MeV) 


Qo (MeV) 


p(n,7)d 


-0.018 


2.224 


rf(p,7)'He 


-0.70 


5.4940 


d{d,p)t 


-0.004 


4.033 


d{d,nfne 
t(d,n)*He 


-0.682 


3.270 


-0.81 


17.590 


^He(d,p)*He 
t(Q,7)"Li 


-0.81 


18.353 


-0.844 


2.467 


='He(a,7)''Be 


-1.74 


1.587 


^Re{n,pf}l 
''Li(p, Q)''He 


0.678 


0.7637 


0.016 


17.3459 


^Be(n,p)^Li 


1.574 


1.6443 



yields change at Vlsh? > 0.025. 

We adopt electromagnetic contributions Mem to nu- 
clear masses based on quantum Monte Carlo calcula- 
tions [2^ , where all electromagnetic contributions 
(including magnetic-moment interactions) are computed 
explicitly. We then multiply Mem by (1 + 5) to vary 
a, and adjust the Q-values accordingly to obtain Q — 
Qo + AMemS- The adopted values of AMem are given 



in Table IB 



In summary, to calculate BBN yields for an altered fine 
structure constant, first apply the changes described in 
BIR. Then apply these additional corrections: 



• Coulomb normalization factor: Multiply 
charged-particle rates by a factor of (1 + 5). 



all 



• Radiative captures: Multiply rates for radiative 
captures by another factor of (1 -I- 15). Multiply the 
rates for '^He(Q!, 7)^Be and ■^H(Q!,7)^Li by the ad- 
ditional factors / in Table ||. 

• Final-state Coulomb interactions: Multiply rates 
with two charged particles in the final state by a 
factor {1 + aS), using the values of a given for each 
reaction in Table |. 

• Masses: Adjust the reaction Q-values using the 
electromagnetic contributions shown in Table III: 
Q = Qo + AMemS. 

Most of these effects have the opposite sign from the 
exponential factor considered in BIR. We consequently 
expect a weaker dependence on a for the reaction rates 
and resulting BBN yields. 
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The uncertainty in this improved treatment is difh- 
cult to assess. The adjustment of the Q-values is rela- 
tively model-independent. Adjustment of the rates, on 
the other hand, is less certain because we have used ap- 
proximations to the Coulomb penetrabilities, and also be- 
cause we have neglected everything happening inside the 
nucleus. The result is certainly better (and closer to self- 
consistency) than the overestimate that results from as- 
suming that ^-factors are independent of a. Our approx- 
imations to barrier penetrabilities almost always overes- 
timate the effect of varying a. The overestimates are 
generally 5-10% of the variation, except where ZxZi = 1; 
here, they may be underestimates by up to 10% of the 
variation. These percentages are relative to barrier pen- 
etration factors for i? = 2 fm and I < 2, and not relative 
to the (unknown and perhaps unknowable) total change 
in the cross section when a changes. The most extreme 
limit to the validity of the method is where states become 
unbound, i.c, Q — Qq + AMem^ = 0. This occurs at 
about S — 0.5 for the first excited state of ''Be, although 
our radiative capture calculation fails for numerical rea- 
sons beyond about S = 0.3. Our treatment of final-state 
barrier penetration also assumes that S is "small" (being 
based on first derivatives), but these corrections turn out 
to be relatively unimportant. We assume symmetric lim- 
its to the validity of our analysis, —0.3 < S < 0.3, with 
less-questionable validity at —0.1 < (5 < 0.1. 



III. LABORATORY RATES 

We note finally that BIR used reaction rates that were 
already out of date, and they did not cite sources for the 
rates they did use. All of their important rates come 
from the evaluation by Smith, Kawano, and Malaney 
(SKM) (and the others trace to the Caughlan and 
Fowler compilations [^). In the time since SKM (1993), 
some of the most uncertain cross sections have been re- 
measured. Newer, better rates exist for the processes 
d(p,7)3He g ^(a'7)^Li (2|, and ^}ie{n,pfU 
The revised rate for the first reaction has the usual func- 
tional form derived by the S'-factor formalism (though 
its parameters are actually fitted to a numerical integra- 
tion), and the expression from BIR, as corrected in Sec. 
II, may plausibly provide its a dependence. We have thus 
substituted the coefficients from Ref. |^ into the BIR 
expression. 

For the second reaction, the revised rate is again a fit 
to a numerical integration, but the functional form used 
does not show up in the BIR calculation. Since we have 
a physical model for the variation of this rate with a, 
and the BIR formalism does in fact account for much of 
the a dependence in this model, we adopt the following 
somewhat convoluted procedure: we first computed a fit- 
ting formula for the ratio of the revised rate to the SKM 
rate. In our BBN code, we multiply the BIR a-dependent 
rate by this fitting formula. We then apply the correc- 
tions from Sec. II. The result provides both the correct, 



up-to-date rate at a = ao, and the a dependence com- 
puted from our direct capture model. The SKM S'-factor 
contains a term of the form 



S{E) = S{0)q{E)e-^^, 



(6) 



which we drop in this procedure to more accurately repro- 
duce the a dependence computed in our model. This is 
not a problem with the SKM rate evaluation, but rather 
a question of producing simple, adequate fitting functions 
for the a dependence of the rate. (Note that in principle, 
this term arises from the external-capture nature of the 
reaction, so that (3 should almost certainly depend on a 
- a possibility not considered in BIR.) 

Finally, we updated the rate for ^He(n,p)'^II, which 
does not need the BIR treatment because of the neu- 
tron in the entrance channel. In summary, our cal- 
culations used the SKM rates, with a variation esti- 
mated as described above, and with updates from Refs. 
p7| , ^ ^ . We compared a calculation using these rates 
and Q!BBN = CKQ with the corresponding results from the 
NACRE [0 and NoUett & Buries Q evaluations of the 
rates (which are not suited to computing a variation). 
The results were in agreement within the uncertainties. 

To compare our calculation with observed abundances, 
it will be necessary to take uncertainties on the labora- 
tory measurements into account. The effects of these 
uncertainties on the computed abundances can, in gen- 
eral, only be estimated by Monte Carlo methods. Since 
we compute a large grid of models in ^Ish^ and a, such 
a procedure would be time-consuming. Instead, we make 
very simple, conservative estimates of the uncertainties. 
The uncertainties on the computed helium mass fraction 
Yp have been estimated by Lopez & Turner to be 
±0.0004 from experiment and ±0.0002 from theory. We 
have not included their elaborate treatment of small ef- 
fects in our calculation, and we will assume an error of 
0.0005 on the calculation. For Qsh'^ > 0.011, the uncer- 
tainty in D/H arises mainly from the rate for d(j>, 7)"^He; 
we use an estimate of 9% here, based on the uncertainty 
in D/H at ^ph^ = 0.036, its maximum in Ref. |3l). Sim- 
ilarly, the uncertainty in ''Li/H at ViBh? > 0.011 is gov- 
erned by the uncertainty in the rate for '^He(Q!, 7)^Be. We 
take as the uncertainty in ^Li/H its maximum of 0.05 dex 
at (S = from Ref. ||l|- We will see below that ^He/H 
depends only weakly on aBBN- Moreover, the identifi- 
cation of the observed "^He/H with the primordial value 
p3| seems to us problematic, so we will not use '^He in 
our examination of the evidence for varying a. 



IV. RESULTS 

We incorporated all of the a dependence discussed 
above into the standard BBN code |34, |35j and com- 
puted yields. Figure |l] shows the resulting light-element 
abundances for S = 0, ±0.05, compared with a calcula- 
tion applying the BIR formalism to updated laboratory 
rates. For D/H and "^Li/H, we find that the a dependence 
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FIG. 1: Light element abundances for 5 = (solid lines), 
S = ±0.05 from the BIR method (dash-dotted lines) and S — 
±0.05 from our calculation (dashed lines), ub and axe the 
number densities of baryons and photons, respectively. 



is weaker than that found by BIR, though the effect still 
works in the same direction. This can be understood 
as follows: the strongest a dependence in their calcula- 
tion is in the suppression of charged-particle rates by the 
negative exponential in those rates. In our calculation, 
positive powers of a from current coupling and careful 
treatment of the barrier penetration partially cancel this 
suppression. The ^Li also receives a somewhat smaller 
contribution from adjustments in outgoing photon en- 
ergy which we modelled. Effects of the changing Q values 
on the reverse reactions become important for ^Li/H at 
large ^Ish^, making our estimate coincidentally equal to 
that of BIR at Q.sh'^ > 0.065. Our estimate of final-state 
effects makes a much smaller contribution to the yields 
than any of the other corrections we have added to the 
calculation. 

Our results for the "'He abundance agree with BIR. 
This abundance depends on a primarily through the 
neutron-proton mass difference ?7i„ — mp which governs 
the neutron/proton ratio at the start of BBN. Like BIR, 
we use the phenomenological prescription of Gasser and 
Leutwyler p6| to model the a dependence of the nucleon 
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FIG. 2; Fractional change in yields from their value at 5 = 0, 
computed at ns/n^ = 5.6 x 10~^° {^Ish^ = 0.0204). Solid 
curves are from our calculation. Dashed curves indicate re- 
sults from the methods of BIR. 

mass splitting: 

TO„ -TOp = 1.3(l + c(5) MeV, (7) 

where c = —0.59. As pointed out in BIR, this model is 
plausible, but certainly not definitive, so the ^He yields 
are less certain than those of other light nuclides. We 
use the Gasser & Leutwyler value because we did not 
find competing, explicit values of c in the literature. The 
abundances of the other nuclides are for the most part 
affected only weakly by this uncertainty, but the nucleon 
mass splitting makes an important contribution to the a 
dependence of ''Li/H at low fisft,^ < 0.010. We show the 
dependences of D/H, ■^He/H, and ^Li/H on 5 at ^Ish^ — 
0.0204 in Fig. |. 

V. OBSERVATIONAL CONSTRAINTS 

Given observationally-inferred primordial abundances, 
our calculations can be used to constrain possible combi- 
nations offish^ and S. In principle, two abundances suf- 
fice to fix values of the two parameters, and a third abun- 
dance provides a check on the result. In reality, there are 
systematic errors associated with both the measurements 
and the varying-a calculation, and we find no solution 
which satisfies the D, ^He, and ''Li observations simulta- 
neously. ('^He/H shows relatively weak dependence on a; 
it may in principle help constrain Vlsh? so that a may 
be extracted from other abundances, given the primor- 
dial ratio.) We now examine the constraints provided by 
the present observations. This discussion is in large part 
supplementary to those in BIR and in Avelino et al. ]37|] 
(who used the BIR rates, but also applied statistical er- 
ror estimation to a subset of the data). However, we also 
draw upon our improved estimation of the yields. The 
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reader should keep in mind that the primordial light- 
nuclide abundances require difhcult observational work 
and are by no means settled questions; we will attempt 
to point out the most serious problems, but the reader is 
referred to recent reviews of the topic for critical discus- 
sion of the uncertainties and appropriate caveats ^ . 
It what follows, we assume Gaussian distributions for all 
uncertainties and include the roughly-estimated yield un- 
certainties of Sec. 



[II. 



A value for a single light-nuclide abundance specifies 
one constraint relating Vlsh? and 5, and therefore defines 
a curve in the {VLBh^,5) plane. As a first example, we 
show in Fig. ^ the constraint in the (fis/i^, i5) plane pro- 
vided by an average of the extragalactic D/H measure- 
ments of Tytler and collaborators, (3.0±0.4) x 10^^ [|o). 
(There exist measurements by other groups, which we ex- 
clude because of concerns about signal-to-noise Q and 
complex velocity structure |4|] ; these tend to increase 
the scatter.) The largest weakness of D/H as a probe of 
BBN is the difficulty of obtaining large enough numbers 
of high-quality data to be confident in their analysis. Be- 
cause there is uncertainty in the observations, there is 
uncertainty in specifying the curve; this is indicated by 
shaded regions delimiting 1, 2, and Scr confidence levels. 
We note that the error estimate of Ref. |Q is rather 
conservative. 

The analogous curves based on halo-star observations 
of ^Li are shown in Fig. ^. We derive these curves from 
the primordial ^Li abundance advocated by Ryan Q| 
and Ryan et al. [||, logio(^Li/H) = -9.9lt|^:i^. It 
should be kept in mind that there could be factor-of- 
two effects associated with Li destruction in the observed 
stars Q and with the stellar atmosphere mod- 

els used to derive abundances ji^. There are two curves 
for ^Li/H, reflecting the two paths for its production in 
BBN: a \ow-^b case in which ''Li is produced directly, 
and a high-f^s case in which ''Li is produced as ''Be, 
followed by beta decay. Again, confidence intervals are 
indicated by the shaded regions. Comparison of Figs. 
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FIG. 5: Confidence levels in the (fis/i^, 5) plane for D/H and 
^Li/H combined, ignoring *He. The contours reflect la, 2a, 
and 3cr confidence levels. The solution for these two parame- 
ters given the two constraints is Q,Bh? = 0.034, 5 — 0.23, but 
the result fails a check against a third constraint. 



and H shows one point where the D /H and ^Li results can 
be reconciled. This point occurs at Q.Bh? — 0.034 and 
5 = 0.23. Although we consider the range of validity of 
our treatment to be —0.3 < ^ < 0.3, there does appear to 
be a second solution at lower flsh'^ and an extreme neg- 
ative S. Note that the best-fit model for D and ^Li with 
(5 = is rejected at 98% confidence. The joint constraint 
is shown in Fig. |^. 

How seriously should this result be taken? A solu- 
tion can always be found for two abundances [assuming 
that the curves in the (OB/i^,a) plane intersect], so its 
validity must be checked with a third abundance. The 
other abundance at our disposal is that of ''He, as ob- 
served in extragalactic Hii regions. Here, we confront 
unavoidable systematic problems associated with the ex- 
traction of abundances from spectral-line data. The two 
largest compilations of data yield discrepant values for 
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0.05 



0.05 



FIG. 6: Regions of la, 2a, and 3a confidence, based on the 
Yp value from OSS Q (top) and IT ^ (bottom). 



the primordial ^He mass fraction: Olive, Skillman, and 
Steigman (OSS) find Yp = 0.234±0.002 (50|Lwhile Izotov 
and Thuan (IT) find Yp = 0.244 ± 0.002^. The two 
resulting curves are shown in Fig. |^. It seems clear from 
these curves that if agreement could be reached concern- 
ing the value of Yp , quite stringent constraints could be 
placed on S irrespective of O^. However, this impression 
is illusory because such constraints involve the size of the 
electromagnetic contribution to the nucleon mass split- 
ting. As discussed above, this number is difficult to cal- 
culate, and subject to large uncertainty. The confidence 
intervals shown include allowances for a 30% Gaussian 
uncertainty in this parameter (the source of the asym- 
metric intervals). 

Comparing Figs. || and ^, we see that the D+''Li so- 
lution above fails the consistency check with ^He. Mak- 
ing this solution work would require either Yp=0.354 - 
substantially larger than even the solar value of 0.275 
- or else c 0.045 (at Yp ^ 0.244) or c = 0.10 (at 
Yp = 0.234), compared to the current best guess of 
c — —0.59. While c is very uncertain, it cannot be posi- 
tive if the electromagnetic force contributes to the baryon 
masses only through self-energies and interactions among 
constituent quarks [p6[ . Put another way, a fit includ- 



ing all three abundances gives a minimum far from the 
D /H -I- ^Li solution (because of the small error bars on 
Yp), but the likelihood at this minimum is low: a fit us- 
ing the IT Yp is rejected at 97.5% confidence (actually 
worse than if (5 = is assumed, since the extra parameter 
removes a degree of freedom), while a fit using the OSS 
Yp is rejected at 98.0% confidence. The exact numerical 
values of these confidence limits are probably not mean- 
ingful, since the Yp error bars in particular are smaller 
than the systematic uncertainties, but the indication is 
in any case that varying a does not describe the data 
significantly better than leaving it fixed. 

There is thus no solution favored by all three observed 
primordial abundances, even at S = 0. Perhaps one of 
the pairwise combinations of abundances provides an ac- 
ceptable solution that requires only plausible systematic 
problems with the third observed abundance. Because 
the D-t-''Li solution requires such a stretch to match the 
^He observations, we can probably rule it out. The 
D+'*He case is interesting because it does not require 
an "extreme" value for S. This is in the sense that 
existing claims for time variation in a based on metal 
lines in quasar absorption systems are on the order of 
10"^ The best fits for BBN with D and ^He only are 
at Qsh^ = 0.020 ± 0.005, S = -Omitomr (IT) and 
nph^ - OmOtomi S = -0-031°;°;^ (OSS; all errors 2a). 
These solutions have the advantage of agreeing with the 
value of riBh^ derived from observations of the cosmic 
microwave background (CMB). However, ^Li/H is dis- 
crepant by a factor of 2-3, so this is no better than the 
standard asBN = Q^o case in its need for a greatly under- 
estimated primordial ''Li/H. 

The final pair of abundances available is ''He-l-^Li. Any 
■*He constraint crosses both curves defined by the ^Li 
constraint, and either of the two "^He data sets could 
be correct. There are thus four different solutions we 
could consider, but it should suffice to say that they all 
favor \5\ < 0.02 and nph"^ < 0.017. The solution least 
discrepant from the D/H data requires D/H=7.3 x 10~^, 
more than twice the current best observational average 
and greater than the highest value presently claimed. 

Thus, given the present observational data for light- 
element abundances, allowing a to vary does not allow 
BBN to fit the data in a way that is convincingly bet- 
ter than in the simplest q;bbn ~ Q^o case. The existing 
discrepancies between the light-nuclide abundances and 
the detailed predictions of BBN are better-explained in 
terms of systematic effects than in terms of varying a. 
The above analysis produces two plausible constraints on 
«bbn: one based on deuterium that derives its strength 
from the simplicity of the D/H measurement (Fig. y), 
and one based on ^He that derives its strength from the 
presumptive strong dependence of Yp on aBBN (Fig. ||). 
The constraint based on '*He which was derived in BIR 
remains the best that one can do to constrain obbn h> 
dependent of any other parameter, with the caveat that 
the nucleon mass splitting is not well-understood, and so 
the constraint could be much weaker than it appears. 
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The constraint based on D/H is perhaps more secure 
because the varying-a D/H yields depend on relatively 
well- understood physics and because the D/H measure- 
ment is in principle relatively free of systematic pitfalls: 
it involves demonstrably pristine material, and it does 
not depend on radiative transfer or ionization correc- 
tions. The deuterium constraint may be expressed as 



-1.68 ± 0.07 at 2ct, but it has the 



disadvantage that one can obtain large variations in a 
merely by adjusting fig. Given the large regions of high 
confidence level in Fig. ^, ^Li does not provide a useful 
constraint by itself. 

We note that the joint constraint from D/H and the IT 
value of Yp was also examined in Ref . |^ . These authors 
obtained a result S — —0.007 ± 0.009 (compared with 
our 6 = — 0.007to!Qi7 for this case), and identical hmits 
on flsh^ to those found here. One would expect differ- 
ing values of d between their analysis and ours: in our 
calculation, D/H depends more weakly on 6, and the up- 
dated laboratory cross section for the crucial d{p, 7)'^He 
cross section is lower than what they used. It appears 
that these two effects have cancelled by chance at the 
best-fit solution. ("^He should behave the same way in 
both calculations, because we have used the same value 
of c). The differing error limits result from a combina- 
tion of the weaker dependence of d{p, 7)'^He on S in our 
treatment, the inclusion of a 9% cross section uncertainty 
on our predicted D/H, and the 30% uncertainty that we 
have assigned to c. 

Finally, an independent constraint on flsh^ now exists 
from measurements of CMB anisotropics. The current 
results are fleh^ = 0.022 ± 0.003 from BOOMERANG 
H, neh^ = 0.022 from DASI ||, and J^s/i^ = 

0.033±0.007 from MAXIMA g] (all la errors). A single 
direct measurement of ^Ish? is a vertical line in Figs. ^ 
^. Values of aeBN based on these measurements and the 
individual light-nuclide abundances are shown in Table 
IV. A firm value of 5 may not be derived from this table, 
because the relative reliabilities of the inferred primordial 
abundances are not known. However, a general trend 
is clear: only the combination "^Li/H -I- MAXIMA or 
the combination of OSS ^Hc with any of the CMB data 
provides any strong indication of (5 ^ 0. For all other 
data, the 2ct limits either include (5 = or nearly include 
it (missing by 10% or less of the best-fit value). 

In principle, the estimation of Ush? from CMB mea- 
surements also depends on the value of a at matter- 
radiation decouphng |55, |5^, acMB- This may be 
included among the parameters to be determined jointly 
from the CMB, or some value may be assumed as a 
prior. The results quoted above and used to derive Ta- 
ble IV are based on the assumption that acMB — c^o- 



The values of acMB and cibbn may have nothing to do 
with each other: a may evolve between these events. 
However, barring entropy increases after BBN, the same 
51s ft.^ today should be inferred from both BBN and the 
CMB. After marginalizing over a, Avelino et al. ||3^ 
find from a combination of the BOOMERANG, DASI, 



TABLE IV: Constraints obtained by combinining individual 
light-element abundances with Qsh^ as inferred from three 
sets of data on acoustic oscillations in the CMB, and from 
the analysis of Avelino et al. |^^, which marginalized over 
a. Numbers shown are the best-fit 5 and the 5 at the ±2(t 
boundaries. 



Abundance 


CMB experiment 


-2a 


6 


+ 2(7 


"He (Izotov) 


MAXIMA 


-0.053 


-0.017 


-0.0028 




DASI 


-0.031 


-0.0083 


0.0045 




BOOMERANG 


-0.030 


-0.0083 


0.0027 




Marginalized 


-0.011 


-0.0052 


0.0054 


''He (OSS) 


MAXIMA 


-0.12 


-0.039 


-0.020 




DASI 


-0.087 


-0.030 


-0.015 




BOOMERANG 


-0.086 


—0.030 


—0.016 




Marginalized 


-0.078 


-0.027 


-0.014 


D/H 


MAXIMA 


-0.051 


0.22 


0.3 + " 




DASI 


-0.18 


0.027 


0.20 




BOOMERANG 


-0.14 


0.027 


0.18 




Marginalized 


-0.17 


-0.038 


0.15 


^Li/H 


MAXIMA 


0.094 


0.22 


0.3 + " 


DASI 


-0.003 


0.18 


0.3 + " 




BOOMERANG 


0.022 


0.18 


0.3 + " 




Marginalized 


-0.011 


0.149 


0.3 + " 



"2cr upper boundary exceeds validity of model. 



0.019 



-0.004 
-0.002 



at la. For 



and COBE data that 

each BBN nuclide, Table DVl shows the effects of combin- 
ing the BBN constraint with this value of fis/i^. The 
results are not markedly different than for the DASI and 
BOOMERANG results independently. This is not sur- 
prising, because Avelino et al. find acMB = c^Oi 
within the uncertainties. 



VI. CONCLUSIONS 

We have improved the calculation of primordial 
abundances in varying-a scenarios by including correc- 
tions previously neglected: normalization of initial-state 
Coulomb penetrabilities, final-state Coulomb interac- 
tions, photon coupling to nuclear currents, photon en- 
ergy and barrier penetration in external direct captures, 
and electromagnetic contributions to the nuclear masses. 
These affect the interactions between composite nuclei, 
so they affect all of the BBN abundances except that 
of ^He. In sum, these corrections can amount to half 
of the effect previously computed for ^Li/H and about 
1/3 of the previous effect on D/H. The most important 
corrections are the initial-state penetrability and current- 
coupling contributions. The nuclear mass modifications 
become important only for Q.Bh? > 0.036. BBN is less 
sensitive to a than previous calculations indicated, but 
the sign of the dependence remains the same. 

Applying our calculation to the data, the most secure 
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constraint we arrive at is logj^Q ^Ish^ — 5 = —1.68 ± 0.07 
{2a), based on the D/H measurements of Tytler and col- 
laborators |M. Combining this constraint with a con- 
straint nsh^^ 0.022 ± 0.003 from CMB anisotropics 
yields 5 — 0.03 ± 0.07 (Icr). A separate constraint, sub- 
ject to both theoretical and astronomical uncertainties 
which are difhcult to assess, is provided by ^He measure- 
ments: the BIR result of approximately \6\ < 0.1 stands 
as the best that can be done here, about the same as the 
D/H + CMB constraint. 

Other recent work on time-varying a suggests that 
between redshifts 1 and 3, a is smaller than uq by 
5 = (-0.72 ± 0.18) X 10^5 § |, i- Sensitivity to de- 
viations from zero at this level cannot be reached with 
BBN, for both astronomical and nuclear-physics reasons. 
Combining this with the large difference in redshifts be- 



tween BBN {z ~ 10^) and the observed absorption sys- 
tems (1 < z < 3), it is doubtful that BBN can shed much 
light on this claim without a specific model for variation 
of a. 



Acknowledgments 

We thank Richard Battye and Marc Kamionkowski for 
helpful conversations, and R. B. Wiringa for providing 
the EM contributions to the light-nuclide masses. REL 
was funded by a PPARC Advanced Fellowship. KMN 
was supported by DoE grant DE-FG03-92-ER40701 and 
NASA grant NAG5-11725. 



[1] J. K. Webb, V. V. Flambaum, C. W. Churchill, M. J. 
Drinkwater, and J. D. Barrow, Phys. Rev. Lett. 82, 884 
(1999). 

[2] M. T. Murphy, J. K. Webb, V. V. Flambaum, V. A. 
Dzuba, C. W. Churchill, J. X. Prochaska, J. D. Barrow, 
and A. M. Wolfe, Mon. Not. Roy. Astron. Soc. 327, 1208 
(2001). 

[3] M. T. Murphy, J. K. Webb, V. V. Flambaum, J. X. 

Prochaska, and A. M. Wolfe, Mon. Not. Roy. Astron. 

Soc. 327, 1237 (2001). 
[4] J. K. Webb, M. T. Murphy, V. V. Flambaum, V. A. 

Dzuba, J. D. Barrow, C. W. Churchill, J. X. Prochaska, 

and A. M. Wolfe, Phys. Rev. Lett. 87, 091301 (2001). 
[5] J. F. Donoghue, hep-ph/0101130 . 

[6] T. Banks, M. Dine, and M. Douglas, Phys. Rev. Lett. 

88, 131301 (2002). 
[7] J. D. Bekenstein, Phys. Rev. D 25, 1527 (1982). 
[8] J. D. Barrow, H. B. Sandvik, and J. Maguijo, Phys. Rev. 

D 65, 063504 (2002). 
[9] E. Kolb, M. Perry, and T. Walker, Phys. Rev. D 33, 869 

(1986). 

[10] J. D. Barrow, Phys. Rev. D 35, 1805 (1987). 
[11] B. A. CampbeU and K. A. Olive, Phys. Lett B 345, 429 
(1995). 

[12] L. Bergstrom, S. Iguri, and H. Rubinstein, Phys. Rev. D 

60, 045005 (1999), abbreviated BIR. 
[13] P. Langacker, G. Segre, and M. J. Strassler, Phys. Lett. 

B 528, 121 (2002). 
[14] T. Dent and M. Fairbairn, astro-ph/0112279 . 
[15] S. R. Beane and M. J. Savage, hep-ph/0206113 . 
[16] X. C. Calmet and H. Fritzsch, hep-ph/0112110 . 
[17] M. Abramowitz and L A. Stegun, eds.. Handbook of 

Mathematical Functions (Dover, New York, 1965). 
[18] J. M. Blatt and V. F. Weisskopf, Theoretical Nuclear 

Physics (Springer, New York, 1979). 
[19] R. F. Christy and I. Duck, Nucl. Phys. 24, 89 (1961). 
[20] K. M. Nollett, Phys. Rev. C 63, 054002 (2001). 
[21] B. Buck, R. A. Baldock, and J. A. Rubio, J. Phys. G 11, 

Lll (1985). 

[22] B. S. Pudliner, V. R. Pandharipande, J. Carlson, S. C. 

Pieper, and R. B. Wiringa, Phys. Rev. C 56, 1720 (1997). 
[23] R. B. Wiringa, private communication. 



[24] M. S. Smith, L. H. Kawano, and R. A. Malaney, Astro- 
phys. J. Suppl. 85, 219 (1993), abbreviated SKM. 

[25] G. R. Caughlan and W. A. Fowler, At. Data Nucl. Data 
Tables 40, 283 (1988). 

[26] G. J. Schmid, B. J. Rice, R. M. Chasteler, M. A. Godwin, 
G. C. Kiang, L. L. Kiang, C. M. Laymon, R. M. Prior, 
D. R. Tilley, and H. R. Weller, Phys. Rev. C 56, 2565 
(1997). 

[27] L. Ma, H. J. Karwowski, C. R. Prune, Z. Ayer, T. C. 
Black, J. C. Blackmon, E. J. Ludwig, M. Viviani, 
A. Kievsky, and R. Schiavilla, Phys. Rev. C 55, 588 
(1997). 

[28] C. R. Prune, R. W. Kavanagh, and C. Rolfs, Phys. Rev. 

C 50, 2205 (1994). 
[29] C. Prune, K. L Hahn, R. W. Kavanagh, and P. R. Wrean, 

Phys. Rev. C 60, 015801 (1999). 
[30] C. Angulo et al, Nucl. Phys. A656, 3 (1999). 
[31] K. M. Nollett and S. Buries, Phys. Rev. D 61, 123505 

(2000) . 

[32] R. Lopez and M. S. Turner, Phys. Rev. D 59, 103502 
(1999). 

[33] T. M. Bania, R. T. Rood, and D. S. Balser, Nature 415, 
54 (2002). 

[34] L. H. Kawano, Fermilab-pub-88/34-a (1988), (unpub- 
lished) . 

[35] L. H. Kawano, Ferrmlab-pub-92/04-a (1992), (unpub- 
lished) . 

[36] J. Gasser and H. Leutwyler, Phys. Rep. 87, 77 (1982). 

[37] P. P. Avelino, S. Esposito, G. Mangano, C. J. A. P. Mar- 
tins, A. Melchiorri, G. Miele, O. Pisanti, G. Rocha, and 
P. T. P. Viana, Phys. Rev. D 64, 102505 (2001). 

[38] K. A. Olive, G. Steigman, and T. P. Walker, Phys. Rept. 
389, 333 (2000). 

[39] D. Tytler, J. M. O'Meara, N. Suzuki, and D. Lubin, 
Physica Scripta 85, 12 (2000). 

[40] J. M. O'Meara, D. Tytler, D. Kirkman, N. Suzuki, J. X. 
Prochaska, D. Lubin, and A. M. Wolfe, Astrophys. J. 
552, 718 (2001). 

[41] M. Pettini and D. V. Bowen, Astrophys. J. 560, 41 

(2001) . 

[42] S. D'Odorico, M. Dessauges-Zavadsky, and P. Molaro, 
Astron. Astrophys. 368, L21 (2001). 



10 



[43] S. A. Levshakov, M. Dessauges-Zavadsky, and S. D. 
D'Odorico, Astrophys. J. 565, 696 (2002). 

[44] S. Ryan, Proc. lAU Symposium 198, "The Light Ele- 
ments and Their Evolution" (Astronomical Society of the 
Pacific, 2000), p. 249. 

[45] S. G. Ryan, T. C. Beers, K. A. Olive, B. D. Fields, and 
J. E. Norris, Astropliys. J. 530, L57 (2000). 

[46] S. Vauclair and C. Cliarbonnel, Astron. Astrophys. 295, 
715 (1995). 

[47] M. H. Pinsonneauh, T. P. Walker, G. Steigman, and 
V. K. Narayanan, Astrophys. J. 527, 180 (1999). 

[48] M. Safaris and A. Weiss, Astron. Astrophys. 376, 955 
(2001). 

[49] M. Asplund, in L. da Silva, M. Spite, and J. R. 
de Medeiros, eds.. The Light Elements and Their Evolu- 
tion (Astronomical Society of the Pacific, Provo, 2000), 



fAU Symp. 198. 
[50] K. Olive, G. Steigman, and E. Skillman, Astrophys. J. 

483, 788 (1998). 
[51] Y. Izotov and T. Thuan, Astrophys. J. 500, 188 (1998). 
[52] C. B. Netterfield et ai, Astrophys. J. (2002), in press; 

astro-p h/0104460| . 
[53] C. Pryke, N. W. Halverson, E. M. Leitch, J. Kovac, J. E. 

Carlstrom, W. L. Holzapfel, and M. Dragovan, Astro- 
phys. J. 568, 46 (2002). 
[54] R. Stompor et al, Astrophys. J. 561, L7 (2001). 
[55] M. Kaplinghat, R. Scherrer, and M. Turner, Phys. Rev. 

D 60, 023516 (1999). 
[56] S. Hannestad, Phys. Rev. D 60, 023515 (1999). 
[57] R. A. Battye, R. Crittenden, and J. WeUer, Phys. Rev. 

D 63, 043505 (2001). 



